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ABSTRACT 

Aims. The knowledge of the properties of the youngest radio sources is very important in order to trace the earliest phase of the 
evolution of the radio emission. RXJ1459-I-3337, with its high turnover frequency (~ 25 GHz) provides a unique opportunity 
to study this class of extreme objects. 

Methods. High-sensitivity multi-frequency VLA observations have been carried out to measure the flux-density with high 
accuracy, while multi-frequency VLBA observations were performed, aimed at determining the pc-scale structure. Archival 
ROSAT data have been used to infer the X-ray luminosity. 

Results. The comparison between our new VLA data and those available in the literature shows a steady increment of 
the flux-density in the optically-thick part of the spectrum and a decrement of the turnover frequency. In the optically- 
thin regime, the source flux density has already started to decrease. Such a variability can be explained in terms of 
an adiabatically-expanding homogeneous radio component. The frequency range spanned by our VLBA observations, 
together with the resolution achieved, allows us to determine the source size and the turnover frequency, and then to 
derive the magnetic fleld directly from these observable quantities. The value obtained in this way is in good agreement with 
that computed assuming equipartition condition. A similar value is also obtained by comparing the radio and X-ray luminosities. 

Key words, galaxies: active - galaxies; evolution - radio continuum: general - magnetic fields - radiation mechanisms: non- 
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1. Introduction 

The radio emission of extragalactic sources is synchrotron 
radiation produced by relativistic electrons with a power- 
law energy distribution. The relativistic electrons produce 
a power-law radio spectrum which reaches its maximum 
in correspondence to the turnover frequency vt- 
At frequencies below vt, the spectrum turns over likely 
due to Synchrotron Self- Absorption (Snellen et al. I2000p . 
although free-free absorption may also play a role (Mutoh 
et al. [^00^ Kameno et al. lMTOl) . 

The evolutionary models proposed to interpret the various 
stages of the life-cycle of radio sources relate the typical 
spectral peak of very small and then compact region with 
age. In the evolutionary scenario, the peak of the spectrum 
progressively moves toward lower frequencies as the ra- 
dio source expands/grows. In this framework, GHz-peaked 
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spectrum (GPS) and compact steep-spectrum (CSS) radio 
sources, the former with vt ^ 1 GHz and the latter with 
^ 100 — 500 MHz, are considered to represent early 
stages in the individual radio source evolution. 
The timescale of the spectral evolution in extremely young 
radio sources is very short, of the order of a few tens of 
years fPallacasa I2003p . This implies that the ideal tar- 
gets to investigate how the radio spectrum evolves, and 
that are the main mechanisms at work, must be sought 
among sources with very high turnover frequencies. The 
existence of "extreme" GPS sources with turnover fre- 
quency above 10 GHz, termed "high frequency peakers" 
(HFPs) by Dallacasa et al. (|2000p , is expected from all the 
radio source evolutionary models. Their detection is, how- 
ever, quite difficult since they are short-lived objects evolv- 
ing into GPS and then GSS sources, and to observational 
limitations preventing their selection. For a source to be 
recognized as an extreme GPS object, the characteristic 
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turnover frequency must lie within the frequency range 
sampled by large area surveys. The frequencies of the sur- 
veys currently available are not high enough to allow an 
efficient selection of high-frequency peaking objects, since 
there is insufficient information above ^ 5 GHz. 
The radio source RXJ1459-I-3337 was one of the first ob- 
jects to be recognized as an "extreme" GPS with the 
turnover frequency at about 30 GHz and a peak flux den- 
sity of about 800 mJy as reported by Edge et al. (|1996p . 
This object is associated with a quasar at z = 0.6448. 
It was detected by ROSAT (Brinkmann et al. 12000)) and 
shows an X-ray luminosity Lx ~ 2.8 x lO^'^erg/s. 
Multi-epoch VLA observations indicate strong variabil- 
ity in the optically-thick part of its radio spectrum. In 
particular, the flux-density variability at 5 GHz (i.e. in 
the optically-thick regime) seems to steadily increase. 
Although flux-density variability is quite rare in GPS 
sources (0'Dea [1998p . an increment of the flux density at 
frequencies below vt may be a consequence of the evolution 
of a synchrotron self-absorbed spectrum in an expanding 
component. 

In this paper we present the results of multi-frequency 
VLA and VLB A observations conducted in 2003 and 2005, 
respectively. In order to identify the mechanisms at the ba- 
sis of the spectral evolution we combine the information 
on the physical conditions provided by our new VLA data 
with that from archival data obtained in 1996 and 1999. 
The frequency range sampled by the new VLBA obser- 
vations, together with their high-resolution, allows us to 
compute the magnetic field by means of observable quanti- 
ties only, such as peak frequency, peak flux density and the 
source size, as from the synchrotron theory (Kellermann & 
P auliny-Tot h II 9 8 D see also Section 4.2). The comparison 
between this value with the field strength obtained as- 
suming equipartition provides important information on 
the physical conditions of the radio source. 

Throughout this paper, we assume the following 
cosmology: Hq = 71 kms~^ Mpc~^, Hm = 0.27 and 
fl\ = 0.73, in a flat Universe. 



2. Observations and data reduction 

Our target RXJ1459-(-3337 was observed with the VLA 
on September 12th 2003 at 8 independent frequencies 
(1.4, 1.7, 4.5, 4.9, 8.1, 8.4, 15.0 and 22.2 GHz) during 
the monitoring program of HFP candidates (Dallacasa 
et al. [20001 Tinti et al. [2005)) . the characteristics of the 
observations and the data reduction are described in 
Orienti et al. p007p . In VLA data the r.m.s noise level 
on the image plane is not relevant (always well below 1 
mJy), if compared to the main uncertainty coming from 
the amplitude calibration errors, which are within (Itr) 
3% at 1.4, 1.7, 4.5, 5.0, 8.1 and 8.4 GHz, and 5% at 15 
and 22 GHz, and which are predominant in the case of 
relatively strong radio sources such as RXJ1459-f 3337. 
VLBA observations of the radio source RXJ1459-f3337 
were carried out on April 4th 2005 at 1.7, 5, 8.4, 15, 22 




Fig. 2. Plot of amplitude vs projected baseline length at 
22 GHz. 



and 43 GHz, in full polarization mode with a recording 
band-width of 32 MHz at 128 Mbps, for a total time of 
10 hours. The correlation was performed at the VLBA 
correlator in Socorro and the data reduction was carried 
out with the NRAO AIPS package. After the application 
of system temperature and antenna gain, the amplitudes 
were checked using the data on 4C39.25 (J0927-f 3902). 
The error on the absolute flux density scale is generally 
within 3%-10%, being worse at the highest frequencies. 
The same source 4C39.25 was used to generate the 
bandpass correction at each frequency. However, at 43.2 
GHz a problem on the scan on the calibrator precluded a 
good bandpass calibration. 

Images of RXJ1459-I-3337 at each frequency were pro- 
duced after a number of phase-only self-calibration 
iterations (Fig. [1]) . The source was found to be marginally 
resolved at all frequencies (Fig. The source flux 
densities at each frequency have been measured on both 
VLBA and VLA images in several ways, giving similar 
results, and the values are reported in Table[TJ The VLBA 
flux-densities at each frequency are consistent within 
< 10% with those measured by the VLA, indicating the 
lack of low-surface brightness features on the pc-scale. 



3. Results 

GPS radio sources are usually characterized by the lack 
of any significant flux-density variability and they can 
be considered as the least variable class of extragalac- 
tic radio sources (O'Dea I1998p . with an average varia- 
tion within ~5% (Stanghellini et al. 12005^ . In the case 
of RXJ1459-I-3337, however, a comparison between our 
new simultaneous multi-frequency VLA observations with 
those available from the literature has pinpointed a sub- 
stantial variability in the optically thick part of the spec- 
trum. 

Fig. [3] shows the light-curves of RXJ1459-)-3337 at each 
frequency: from these plots it seems that at 1.4 and 5 
GHz, well below the turnover frequency, the flux density 
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Date 5*1.4 Si. 7 S4.5 S5.0 Ss.i Ssa S\5 S22 S43 

mjy mjy mjy mjy mjy mjy mjy mjy mjy 

VLA Sept 2003 25 31 194 221 410 423 521 356 

VLBA Apr 2005 33 221 446 519 348 198 



Table 1. The VLA and VLBA flux densities of RXJ1459+3337. 



has steadily increased, while in the optically-thin regime 
(i.e. 22 GHz), the flux density has decreased. Following the 
approach from Dallacasa et al. (|2000p . we fit the simulta- 
neous radio spectra with a purely analytical function used 
to determine the peak fiux density Sm and the frequency 
Vm at which it occurs: 

Logs' ^a- y/b^ + (cLog(j/) - d)2 

In this equation, the parameters a, 5, c and d are purely 
numeric, and do not provide any direct physical in- 
formation. Based on the results of the fit, we found 
that the spectral peak has moved to lower frequencies, 
from -24±1 GHz in 1996 to 17±2 GHz in 1999 and to 
12.5±0.5 in 2003 (Fig. |4]). The errors on these quantities 
were calculated following the error propagation theory. 
The decrement of the spectral peak together with the 
steady increment of the flux-density in the optically-thick 
part of spectrum suggest that the source is adiabatically 
expanding. Despite the high resolution achieved by VLBA 
observations, the radio source RXJ1459-I-3337 is only 
marginally resolved even at 43 GHz. Since the inferred 
source size is substantially smaller than the VLBA beam 
size at the highest frequency, our measurements are not 
sensitive to small changes on the angular size and thus 
cannot confirm the source expansion suggested by the 
spectrum evolution. 



4. Discussion 

The strong variability in the optically thick part of the 
spectrum is not common among sources peaking below 5 
GHz (O'Dea [1998p . However, the timescale of evolution 
of extreme GPS sources appears to be quite short, of the 
order of a few tens of years. In the context of the youth 
scenario, such a result is not unexpected since very young 
radio sources should evolve quite rapidly (Dallacasa 
[2003)1 . 

In the following discussion we investigate the physical 
conditions of the radio emission and which mechanisms 
are involved in the evolution of the radio spectrum. 
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Fig. 4. The VLA radio spectrum of RXJ1459-h3337 
in 1996 {squares, Edge et al. I1996p . in 1999 {crosses, 
Dallacasa I2003p . and in 2003 {triangles, these new data). 



4.1. The evolution of the radio spectrum 

As noted above, the radio emission in extragalactic radio 
sources is due to synchrotron radiation from relativistic 
particles with a power-law energy distribution. However, 
in very small objects the observed radio spectrum signifi- 
cantly departs from the power-law shape, turning over at 
frequencies below the peak. Such a deviation is mostly ex- 
plained in terms of synchrotron self-absorption. 
The mechanism which plays the major role in the evo- 
lution of the optically-thick part of the spectrum is the 
adiabatic expansion. In this regime we have that 

S{v)ozH-Kh^ (1) 

(Kellermann & Pauliny-Toth I198ip where S is the flux 
density at a given frequency v in the optically-thick part 
of the spectrum, H is the magnetic fleld and 6 is the 
angular size of the emitting region. As the radio source 
adiabatically expands, the opacity decreases, the turnover 
frequency moves to lower frequencies and the flux density 
increases at frequencies below the turnover frequency. 
Is it possible to relate the flux-density variability observed 
in RXJ1459+3337 to adiabatic expansions? 
First of all, we assume that the radio emission is due 
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to a homogeneous component which is adiabaticaUy 
expanding at a constant rate: 

where Oq is the source angular size at the time to and 9 
at the time to + At. 

We also assume that the magnetic field is frozen in the 
plasma: 

\to + AtJ 

where Ho is the magnetic field at the time to and H the 
magnetic field at the time to+At. From Eq.[l]we find that: 



where vo and v are the turnover fi'cquency at the time to 
and to + At respectively. 

We estimate the epoch of the origin of radio emission 
by fitting the frequency peak measured at the different 
epochs with the same function used to fit the multi-epoch 
flux density. In this case, y=(j/o/i')^/*, x = t ~ to and 
a the age the radio emission had at the epoch to- vo is 
the peak frequency at the oldest epoch available (^o, in 
this case 1996). From the best fit we find a ~ 47 ± 10 
years, implying that the radio emission originated at 
^ 1949 ± 10, still in agreement with the values obtained 
from the fiux density increment, although with higher 
uncertainties, since the fit has been performed on 3 data 
points only. The good agreement on the determination 
of the source age may indicate that the radio source is 
evolving with the physical conditions predicted by simple 
self-similar growth models, as previously assumed. 



s_ _ r tp + At 

So V to 



(2) 



where So and S are the fiux densities measured at the 
same frequency at the time to and to + At respectively. 
In the case of the 5 GHz light-curve, where enough data 
points are available, we estimate the time of onset of 
radio emission by fitting the fiux density measured at the 
different epochs with the function 



y 



(x/a) 



where y = (S/So) 



1/3 



t — to, a the age the radio 



emission had at epoch to and b a free parameter which is 
almost equal to unity in our fits. As So we consider the 
flux density measured at the oldest epoch available (to', i.e. 
1986). From the best fit to the data we obtain a ~ 30 ± 5 
years, i.e. the birth of the radio emission should have oc- 
curred in 1956±5. In the case of the 1.4 GHz light-curve, 
the availability of only 3 data points makes the determi- 
nation of the age less accurate with a ^ 27 ± 15 years, 
i.e. the radio emission should have started in 1966±15, 
in agreement with the 5 GHz determination. We do not 
perform fits on the light-curves at higher frequencies since 
they are affected by contamination due to the flattening 
occurring near the spectral peak. 

As previously mentioned, adiabatic expansion also 
causes the spectral peak to shift at lower frequencies. We 
know that the turnover frequency is 



(3) 



where the energy of the relativistic particles is defined as 
E oct^^ (Pacholczvk fT970|) . This implies that 



to 



At 



to 



4.2. The magnetic field 

Direct measurements of the magnetic field are very 
difficult to carry out. The magnetic field H can be deter- 
mined from observable quantities, such as the turnover 
frequency ut in GHz, the peak flux density S'max in Jy and 
the source angular sizes 6'niaj and ^min in mas as directly 
inferred from the observations 



In this case we have 



H - 10" 



g2 n2 5 2 
-^maj min t '^max 



(1 + ^)- 



(5) 



(Kellermann & Paulinv-Toth I19ST|) where H is in Gauss. 
From the fit to our new VLA data (see Section 2), we 
obtained i^t ^ 12.5±0.5 GHz and S'max ^550±15 mJy. In 
addition, by using the task JMFIT we measured a source 
angular size of 0.36x0.14 mas {9 = 1.8 x 9fit where 0fit is 
the deconvolved full width at half-maximum of a Gaussian 
brightness distribution as used in the fitting procedure) 
on our best VLBA image, namely that at 22 GHz, since 
the longest spacings at 43 GHz are poorly sampled and 
with a rather high noise. The resulting resolution at 43 
GHz is therefore no better than that at 22 GHz. With 
these parameters we obtain H ^ 0.16 ± 0.03 G. 
Another way to constrain the magnetic field is to assume 
that the radio source is in equipartition conditions. From 
the minimum energy conditions (Pacholczyk I1970p . the 
equipartition magnetic field is 



24 

where the minimum energy density Mmin is 



(6) 



= 7-10" 



(l + fc)4/^ 



and L is the luminosity calculated at 22 GHz, i.e. in the 
optically-thin part of the spectrum, and V the volume 
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of the radio source in pc'^. For the calculation of V, 
we assume that the radio emission has an ellipsoidal 
geometry and an average optically thin spectral index of 
0.7 (S oc Furthermore, proton and electron energies 

are assumed to be equal (k=l), with a filling factor of 
unity (i.e. the source volume is fully and homogeneously 
filled by relativistic plasma). From Eq. [6l we calculate 
an equipartition magnetic field Hcq ~ 0.16 G, in very 
good agreement with the value directly derived from the 
observed spectral quantities. 

Since RXJ1459+3337 is also an X-ray source, we can 
try to constrain the magnetic field by comparing the 
radio and X-ray luminosities, if we assume that all the 
X-ray emission is due to Comptonization of the electrons 
responsible for the synchrotron emission in the radio 
band. In this case we have: 



L: 



ssc 



STTWf 



(7) 



(e.g. Singal I1986P where isyn and Lssc are the syn- 
chrotron and the synchrotron self-Compton luminosity 
respectively and Wf is the radiation field, defined as 



Wf 



3Ls 



yn 



Ac 



with A the surface area of the synchrotron emission, and 
c the speed of light. The ROSAT X-ray flux density of 
RXJ1459-h3337 is 1.56xl0-i2 erg s^^ cm^^ (Brinkmann 
et al. I2000p . which corresponds to an X-ray luminosity 
i'SSC ~ 2.8x10''^ erg/s. From the radio fiux density we 
derive a synchrotron luminosity isyn ~ 1.2x10"*^ erg/s. 
If in Eq. [7] we introduce these luminosities and the source 
emitting area A ~ 3.9 pc^ as derived from VLBA images 
assuming an ellipsoidal geometry (see above), we obtain 
a magnetic field H ^ 0.19 G, still in good agreement with 
the values derived by the previous methods. 
These results strongly suggest that the radio emitting 
plasma in such an extreme object is consistent with the 
equipartition conditions. 



5. Conclusions 

We have presented new multi-frequency VLBA and VLA 
observations of the radio source RXJ1459-I-3337. Given 
the position of the spectral peak at very high frequency, 
this object is classified as an extreme GPS radio source, 
also known as HFP. 

By comparing our VLA observations with data from 
the literature spanning about 17 years, we find that the 
flux-density in the optically-thick part of the spectrum 
has been continuously increasing, while the turnover fre- 
quency has been moving toward lower frequencies. Such a 
regular variability, although not common in GPS sources 
peaking below 5 GHz, can be described in the case the 



radio emission originates from an adiabatically-expanding 
homogeneous component. 

If in this context we compare either the flux-densities at 
the same frequency measured at different times or the 
turnover frequencies, we constrain the age of the radio 
emission, which is about ~50±10 years. 
The resolution achieved by the VLBA observations, 
together with the frequency range sampled, allowed us to 
"observe" the turnover frequency and to directly measure 
parameters, such as the peak flux density and angular size 
of the source with great accuracy. In this way we estimate 
the magnetic field to be 0.16 ± 0.03 G by means of 
observable quantities only. The consistency of this value 
with the field derived for a minimum energy condition 
strongly supports the idea that such extreme objects, 
as RXJ1459-I-3337 may already be in equipartition. The 
availability of ROSAT X-ray observations allowed us to 
infer the magnetic field by comparing radio and X-ray 
luminosities. Even in this case, the value inferred for the 
magnetic field agrees with the previous results. 
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Fig. 1. VLBA images of the source J1459+3337 at 1.7, 5, 8.4, 15, 22 and 43 GHz. On each image we give the fohowing 
information on the plot itself: a) Peak flux density in mJy/beam; b) First contour intensity (mJy/beam), which is 

generally equal to 3 r.m.s. of the noise level; contour levels increase by a factor of 2; c) the dynamic range (the ratio 
of the peak brightness and la); d) The restoring beam, plotted on the bottom left corner of each image. 
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Fig. 3. Light-curves of RXJ1459+3337 at 1.4, 5, 8.4, 15 and 22 GHz. In addition to the data points presented in this 
paper we include those concerning the epochs 1986 (87GB, Gregory et al. 1996), 1991 (Neumann et al. I1994p . 1993 
(Laurent-Muehleisen et al. (fT997|) . 1996 (Edge et aL fTMe]) and 1999 (Dallacasa ,2003^. 



